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Shearlet approximation of functions with discontinuous derivatives 

Philipp Petersen 


Abstract 

We demonstrate that shearlet systems yield superior A-term approximation rates compared with 
wavelet systems of functions whose first or higher order derivatives are smooth away from smooth dis¬ 
continuity curves. We will also provide an improved estimate for the decay of shearlet coefficients that 
intersect a discontinuity curve non-tangentially. 


1 Introduction 


In applied harmonic analysis one important field of study is the design of suitable systems (</3n)nGN C L 2 (M^) 
that can efficiently represent functions / S L 2 (R^) in the sense that / = J2neN^n{f)>fn for a sequence 
(c„(/))„gN- One particularly desirable feature of such a system is that by using only few elements it already 
yields decent approximations of functions taken from some subset of L 2 (R^)- The quality of approximation 
within such a system is customarily measured in terms of the error of the best N- term approximation. 

The best N-term approximation of / is given by 


CTJV(/) 


inf 

E]S! CN, I-Eat |— 


ii/-hvii 


2 

L2(R2)- 


/n= V 


If for some class of functions 0 C L 2 (R^) one has aN^f) = 0{g{N)) for TV —)• oo for all / S 0 and some 
function g, we will call g a best N-term approximation rate for 0. 

One particular type system used efficient representation of functions are wavelet systems, see [6] 

and the introduction in Subsection |2.2[ Wavelet systems have established themselves as a standard tool for 
image analysis. The wavelet construction is based on isotropic scaling of a generator function, which implies 
that the supports of all elements <pn have the same aspect ratio. Due to this construction wavelets perform 
poorly when representing functions that contain anisotropic components. For instance a function could have 
a discontinuity along a smooth curve. In this event wavelets yield non-optimal approximation rates as we 
will recall in Subsection 12.31 

To overcome this shortcoming shearlets were introduced in m- These systems constitute an excellent 
tool for the approximation or representation of natural images due to the fact that they provide almost 
optimal best iV-term approximation rates for functions that have discontinuities along smooth curves. We 


will give more details on shearlet systems in Subsection 2.2 


Naturally one can pose the question about approximation rates of wavelets and shearlets for other classes 
of functions which exhibit anisotropic structures. Let, for example, u be the solution of the operator equation 


Lu = /, 

where the data / exhibits a discontinuity along a smooth curve and L is an elliptic differential operator. 
The approximation rates of shearlets for such functions have not been studied yet, despite the fact, that 
good approximation properties for such functions are crucial for the design of optimal adaptive solvers for 
elliptic partial differential equations. In fact, it is well established, see nmiii], that using wavelet systems 
as ansatz functions, certain partial differential equations (PDFs) can be solved with an adaptive strategy 
with a computational complexity depending on the best iV-term approximation rate of the solution of the 
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PDE admitted by the chosen system. In view of applications of shearlet systems for the discretization of 
PDEs it is thus important to know whether the best A^-term approximation rate of shearlets outperforms 
that of wavelets also for classes of functions that contain solutions of elliptic partial differential equations. 
In fact, if this were not the case, it would not be worthwhile to study adaptive frame methods with shearlets 
instead of wavelets. 

Motivated by these considerations we examine functions with first or higher order derivatives which are 
smooth away from smooth discontinuity curves. We will first consider functions which have a /-th order 
derivative which is cartoon-like in the classical sense, i.e., it is twice continuously differentiable apart from a 
twice continuously differentiable discontinuity curve. Cartoon-like functions are a standard model in image 
processing to model natural images. They where first analyzed in [3]. By now analyzing the approximation 
rates for the class of cartoon-like functions was established as the standard benchmark problem to determine 
the quality of a representation system, such as shearlets [Hill], but also curvelets, [T], and more general 
systems, [TO] . 

We will observe, that shearlet s yie ld a superior best A^-term approximation rate of such functions of 
for N ^ oo (Theorem |3.2| ), whe n co mpared to wavelets whose approximation rate cannot be 
faster than for N ^ oo (Theorem |3.l[ ). However, we will see, that not the discontinuity but the 

the regularity assumption on the smooth parts in the cartoon-like model limits the approximation rates. 
Hence we introduce a smoother cartoon-like model in Subsection 13.31 and we will observe in Theorem 13.31 
that the approximation rates by shearlets improve drastically to for N —>■ oo, while those of 

wavelet systems remain unchanged at for N —>■ oo. In order to obtain the approximation rate of 

Theorem |3.3| we prove an improved estimate for shearlet coefficients that do not intersect a discontinuity 
curve tangentially in Lemma |3.4| 

Contemplating the results of this paper in the context of adaptive frame methods discussed above, it is 
fair to say that the fact that shearlets yield improved approximation rates over wavelets for functions with 
cartoon-like derivatives can serve as a justification for the study of adaptive frame methods with shearlet 
systems instead of wavelets systems. We expect that these results will trigger the developement of new 
adaptive frame methods based on shearlet systems. 

2 Preliminaries 

2.1 Notation 

For 0 < p < oo, we will denote by Lp the usual Lebesque spaces of p—times integrable functions and by 
ip we denote the corresponding sequence spaces. For 1 G N we denote by Wj, the spaces of Z—times weakly 
differentiable functions with derivatives in Lp and by the spaces of Z—times continuously differentiable 
functions. For a set H C we use dB to describe its boundary and xb is its characteristic function. For 
two functions h,g we write Zi < p is there exists a constant C such that Zi(x) < Cg{x) for all x in the domain 
of Zi,p. 

2.2 Wavelets and shearlets 

2D Wavelet systems are constructed from integer translations and dilations of a generator function w G 
L 2 (K^), i.e. for j G N and m G 


Wj,m(x) := 2^uj{2^x — m) for all x G 

It can be shown, that there exist G L2(IR^) such that the wavelet system 

W := m G Z^} U {w* j > 0, m G Z^ i = 1, 2, 3} 

yields an orthonormal basis for L 2 (M^), see [6]. We see that the scaling of the wavelet elements is isotropic. 
In the construction of shearlets this isotropic scaling is replaced by parabolic scaling in combination with a 
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shearing matrix. In particular, shearlet systems are constructed using the following two matrices: 

^2^ 0 


Sk = 


1 k 
0 1 


, and Aj = 


0 2! ' ’ 


where £ Z. 


We use cone-adapted shearlet systems which are defined as follows. 

Definition 2.1. fT^ Let '0 G c = [ci,C 2 ]^ G with ci,C 2 > 0. Then the (cone-adapted) shearlet 

system is defined by 

Sn{(j}, 0, 0, C) = $(0, Cl) U ^-(0, c) U ^(0, c), 

where 

$(0, Cl) := {0o,o,m.o = 0(- - Cim) : m G Z^} , 

^'(0, c) := = 2 ti/,(5'^,Aj. .-M^m) : j G No, |fc| < 2r2l ,m G Z^l , 

1^(0, c) := _i = 2 ti/;(5'Ja^-.-M am) : j G No, |fc| < 2r5l ,771 e Z^l , 


with 'lp{xi,X 2 ) = 'lp{x 2 ,Xi), 

AL := 


% M- - 

0 C2<’ 


C2 0 

0 Cl 


, and Aj = 


24 0 

0 23 


For cone-adapted shearlet systems we will employ the index set A := {(j, fc,m, t) : \i\j > j > 0, |fc| < 
|7|24,m G Z^,t = {1,0,—!}}. Let furthermore (i[_,r/ 2 ,ir/ 2 ] denote the following metric on the torus T = 
[-7r/2,7r/2]: 

'^[- 7 r/ 2 , 7 r/ 2 ](a: b) '.= min(|a — 6|, tt -|- minja, b} — maxja, 6}). 

A shearlet is called separable, if 0(a7i,a:2) = 0^(a:i)0^(a:2) for two functions 0^,0^ G L 2 (K)- 

Under certain assumptions, see m cone-adapted shearlet systems can form a frame, i.e. there exist 
0 < Cl < C 2 < c» such that 

C1II/II2 < ^ C2||/||2- 

AgA 

This implies, that there exists a dual frame (0^)agA such that 

T^:f^L2{R^), (ca)agA ^ T"((cA)AeA) := 5] CA0^, (1) 

AgA 

is a bounded operator and / = T‘^(((/, 0 a))aga) for all / G L 2 , see [5]. 


2.3 Approximation of cartoon-like functions 

The improvement of shearlets over wavelets becomes evident, when one considers their approximation rates 
for classes of functions that model natural images, so-called cartoon-like functions. 

Definition 2.2. The set of cartoon-like functions is given by 

■= {/ e L 2 (K^) ■ f = gi +Xb52, where 51,52 G C'^,supp 51,52 C (0,1)^ 
and B C (0,1)^ with dB G and dB has curvature bounded by c}. 

We call dB the discontinuity curve of /. 

It can be shown that wavelets only achieve a best A^-term approximation rate of 0{N~^) for the class 
of cartoon-like functions, [T]. On the other hand in m and m it was shown that cone-adapted shearlet 
systems achieve a best Wterm approximation rate of 0{N~^ log(7V)^), which - up to the log factor - is the 
optimal approximation rate that any system can achieve, see [3]. 
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3 Results 


We are concerned with approximation rates of functions that have some first or higher order derivative which 
is cartoon-like. Let / € N, and 


a 


:= ■\u G and u°‘ S £^(v) for all |a| = I 


Note that in the definition of cartoon-like functions 52 = 0 is possible, so that not all Z—th derivatives of 
u G £^’^(v) need to have a discontinuity curve. We will show, that the best Wterm approximation rate for 
£^’^{v) by wavelet systems is bounded from below by g{N) = and that there exists a best TV-term 

approximation rate by shearlets of the order of g{N) = 


3.1 Wavelet approximation rates 

In order to analyze the approximation properties of wavelets we measure the smoothness of the function to 
be approximated in a Besov scale. We use the following definition of a Besov space, which can be found in 

m- 

For h > 0 and / G L 2 (R^) we define Ahf = fix + h) — fix). For s = k + r, where k gN and 0 < r < 1 
and 1 < p, g < 00 the Besov space Bp ^ is defined by 


B;,, : = {fG 


< 00 } 


ll/ll 


Bf. 



\h\-n.+qr 



1 

q 


Many wavelet systems on a domain on admit a characterization of Besov spaces ^ by non-linear 
approximation rates in the following sense: 


<oo^uG Bg^, 


( 2 ) 


for 1 /g = s /2 -|- 1/2 and r > s, see for instance [Sl|3]. 

Using § we can now find an upper bound for the approximation rate of wavelets for functions in f 

Theorem 3.1. Let I G N, 0 < p < 00 , u G Wp such that for some |q;| = I we have 


dx[ 


-U = XD, 


where D is a bounded subset o/K^ with smooth boundary curve dD that has bounded curvature. Let W be a 
wavelet system such that (H holds. Then for all e > 0 we have that Cf^iu) ^ 0(TV 1^+^) '^). 

Proof. By a simple computation we obtain that xd ^ B^ ^ for any 0 < g < 00. Consequently, u ^ B]^^ 
for any 0 < g < 00. If (TAr(u) < for iV —00 for some e < 0 it would follow with s = T -I- 1, 

1/g = (s -b l)/2 applied to that u G which is a contradiction. □ 


3.2 Shearlet approximation rates 

Now let us prove that for 0 < 1 G N, shearlet systems obtain a higher approximation rate than wavelets for 
functions in 

Theorem 3.2. Let I > 0, and 0 <v < 00 , and u G and let '0) V'j c) = ('0 a)agA be a shearlet 

system with shearlet generator = (g~) & such that 9 G ^ 2 ( 1 ^^) has compact support, is separable and 
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r*/) 1^(01 ^ min(l, | 6 r)niin(l, l^il '^) min(l, |^ 2 | ") and 

(n) 14^(01 <IMCi)(i + |)-"l, 

where t > b,v > 4, h G Li(K). Further assume that forms a frame. Then 

11“ - XI ^ (3) 

X^E]S! 

where Em C A contains the indices of the N largest coefficients \ {u,'tfx)\ and is the canonical dual frame 
element of tpx ■ The log term in @ can he dropped if I > 1. 

Proof. For I = 0 the result is Theorem 1.3 of [Hj. In the sequel we assume Z > 1. Let 71 be the discontinuity 
curve of Ml := ^ 4 ") “ T'-i discontinuity curve of m_i = “• denote for t G % the 

outer normal at t of 7 ^ by n*(t), * = 1, —1. Now we decompose our index set A. First of all we denote by 

A° = {(j, k,m,i) G A : L = 0}. 

By the compact support of u and ip we have that fp{X € A° : (m, ^ 0} < 00 . Furthermore we denote 

:= {(j, fc,TO, t) e A\ A° : supp n |J 7 ^ = 0}. 

.={- 1 , 1 } 

We denote by A^’“ the indices in A^ := A \ (A° U A^), such that for (j, 2l/^s , m, i) G A there exists some t G 
supp V'j, 2 J/ 2 s,m,iLl 7 t such that we have s'" G (tanpo —3-2“l/^,tanpo+3-2“l/^) andn''(t) = ±(cos po) sinpo) for 
some po G (— 7 r/ 2 , 7 r/ 2 ). Additionally, there are non-aligned shearlet elements that intersect the discontinuity 
curve, we will call these indices A^’^. We have A^’^ = A \ (A° U A^ U A^’“). 

Now let us observe the sizes of the index sets for fixed scale j and we denote these index sets by 
Aj, A^'°', A^"^. We have that fp{{j,k,m,b) G Aj : {u,ipj^k,m,i) 7 ^ 0} < 2^'’ . Observe that, due to their 
bounded curvature, 71 and 7 _i have finite length. Hence we observe that 

|A2’“|<24, \Af\<2iF 

Part 1: A°: 

There are only hnitely many indices in A*^, i.e. we certainly have 

n>N 


where c(m)* denotes a non-increasing rearrangement of {\{u,ipj^k,m,L)\'^)(j,k,m,i.)£A°■ 

Part 2: supp ippk,m,L O U.={-i,i} 7. = 0 : 

In this case the scalar products {u,ipj^k,m,i) decay as if the function was in (7*+^. In other words, by 
invoking Proposition 2.1 of [2] (which was only stated for I = 0, but the extension is straightforward) this 
means that 

n>N 

where c(m)* denotes a non-increasing rearrangement of {\{u,ipj^k,m,L)\'^){j,k,m,L)GA'i ■ 

Part 3: supp ippk,m,L O U.={-i,i} 7. 7^ 0 : 

Using partial integration and the compact support of 9 one obtains 


I (P-! Ipj,k,771,1) I 


5 , 


.u. 


7 : ^j,k,7171,7) 


(5) 


Invoking Proposition 2.2. in [14], which are the standard estimates for shearlets i.e. {I = 0), in combination 
with yields that for (j, k, m, l) G A^’“ 
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and for {j, k, m, l) S we have 


Using p 


2/(1 + 3) and the sizes of 


, A^’^ we can compute, that 


A2 


— ^ ' I ('W, V'i.fc.m.i,) 1^ + ^ ^ I ('W,'0j,fe,m,t) 1^ 

A2,a 


<^25'^2"P(j+')'^'+^252"P(j+')'^'<00. 

jGN jGN 


(6) 


Stechkin’s lemma, see e.g. [7], states that 


(7) 

n>N 

if (c^n)n is a monotonically decreasing sequence such that for s = l/p — 1/2. Using Stechkin’s 

lemma, we obtain that X)n>Af''(^)n ^ lV“(*+2)^ where c(rt)* denotes a non-increasing rearrangement of 
(|(u,'!/'i,fc, 7 Ti,t)P)(j,/c,m,t)GA 2 - Using the estimates from Part 1,2, and 3 and invoking the dual frame charac¬ 
terization Q yields 

||u- H c(m)* < 

where c(u)* denotes a non-increasing rearrangement of (|(m,'!/'a)P)agA- □ 



Figure 1: Left: Function u with cartoon-like derivatives. Middle: Cartoon-like first derivative of u in 
vertical direction. Right: Plot on a logarithmic scale of the predicted decay of 0{n~‘^) for n —oo due to 
^ (dashed red line) and actual decay of (c(u )*)2 (solid blue line). For better 


We give a numerical example to illustrate Theorem 3^ We depict in Figure the shearlet coefficients 
of a function that has a cartoon-like derivative in a:: 2 -direction. By these coefficients should decay as 
0{N~‘^) for N ^ oo. In fact the actual decay matches very closely the predicted decay of 0(N~‘^) for 
N ^ oo. The experiments where made with a subsampled version of the shearlet decomposition algorithm 
of ShearLab, m- 


3.3 A different cartoon model 


We saw in the proof of Theorem 3.2 in Q that the main bottleneck for the decay of the shearlet coefficients 
is Q due to the regularity of the functions gi,g 2 of the cartoon-like function / = gi+XDg 2 - We can examine 
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what happens when we assume more regularity. To keep technicalities at a minimum, we only consider the 
case where 51,32 G C°°, but other cases can be studied similarly. We study the following cartoon model: 
Let 1 < Z G N, and 


:= <{uG 


| 2 ) . ^u = ' 
' dx°^ 


and G for all |q!| = I 


where 


^{v) = {/ G L 2 (M^) : / = 31 +XB32, where 31,32 G C°°,supp 31,32 C (0,1)^ 
and B C (0,1)^ with dB G and dB has curvature bounded by u}. 


Obviously C From Theorem 3.1 we know that the best A^-term approximation rate of 


wavelets for this class of functions is bounded from below by 

We can now state the approximation rate of shearlet systems for the class In the proof we will 

make use of Lemma |3.4| which is given subsequent to the following theorem: 

Theorem 3.3. Let I > 1, and 0 < v < oo, and u G £^'°°{v) and let = (^a)agA &e a shearlet 

system with shearlet generator ijj = & such that with L G N, L > Z + 29/4, 0 G C'^(K^) has compact 

support, is separable and 

h) 1^(01 ^ min(l, |^i|^)min(l, |^i|-'^)min(l, 1 ^ 21 "“') and 

(n) + 

where r > 5, b > 4, Zi G Li(IR). Assume that 9 has M > 2L vanishing moments in xi-direction. Further 
assume that (■!/a)agA forms a frame. Then for every e > 0 


u - 


Ae 


( 8 ) 


where Em C A contains the indices of the N largest coefficients |(u,'0a)| and is the canonical dual frame 
element ofipx. 

Before we can present the proof of Theorem |3.3| we require the following improved estimate for shearlet 
elements that intersect a discontinuity curve of a cartoon-like function non-tangentially: 

Lemma 3.4. Let if = ® ff be a separable shearlet where G H ^ 2 ( 1 ^^) are compactly 

supported and has M G N vanishing moments. Let u = gi + Xb 92 with 31,32 G and B C (0,1)^ 

with 7 = dB G and 7 has bounded curvature. Let R> L and L + P < M. 

Let {j, k, m, l) G A such that there exists t G supp H 7 such that the normal n(t) of 'j at t obeys 

n[t) = ±(cos po) sinpo) for some po G [— 7 r/ 2 , 7 r/ 2 ] and cZ[_,n./ 2 , 7 r/ 2 ](arctan(( 2 “-Z/^fc)''), po) > /3 > 0. Then we 
have 

I I < C„(2-(^+1H + 

where Cu is a constant depending only on u. 

Proof. Assume w.l.o.g. that m = 0, k = 0, l = 1, for general |fc| < 2-^/^ one can apply a transformation to 
revert back to Zc = 0. The proof for the case i = — 1 is identical to (. = 1. Since k = 0, l = 1 we have that 
Po ^ (—/3, /3) and hence locally in a neighborhood of t = (ti, ^ 2 ), there exists i? : M —)• K such that 7 is given 
by Xi !->■ (xi -I- ti,E{xi) + < 2 )- E is bounded in a neighborhood of t independent of the chosen t since its 
slope is bounded due to the constraint on po. Furthermore E is i?-times continuously differentiable. Now 
we can apply the transformation theorem to the above equation to obtain 

{'fj,0,0,l,u) = / tpjfip^i{xi,X2- E{xi))u{xi,X 2- E{xi))dx. 

JK2 
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The discontinuity curve of u is locally given by x i—>■ (xi +< 1 ,^ 2 ) and hence u(xi, X 2 — E{xi)) is for fixed 
X 2 - We can also assume, that for some c > 0 we have supp V'j. 0 , 0,1 C A~^{[—c,c\^). Hence we can replace 
(xi,X 2 ) I—>■ m(xi,X 2 — E{xi)) by a P-th order Taylor approximation U to obtain 


/ '^j,o,o,i{xi,X2-E{xi))u{xi,X2-E{xi))dx- / ?/'j_o,o,i(a:i,a;2-P(xi))C/(xi,X2)dx| < 2 


A-\[-c,cY) 


Now we decompose V'j,o,o.i(a^i)— E{xi)) = (x 2 — £’(xi))). Let E[{-,X 2 ) be a L-th 

order Taylor approximation of xi 1 —>■ (jA{2^ {x 2 — E{xi)))), then 


Iaj\[-c,c]^) 


22'0^(2'^xi)2i^^(22 (x2 — E(xi)))U{xi,X2)dx 
2^{2^xi)2iH{xi, X2)U{xi, X2)dx + 


Ja-\[-c,c]^) 

By construction ElU is a polynomial of order L + P. Since 'tp^ has M > L + P vanishing moments 

I 2^ip^{2^xi)2iH{xi, X2)U{xi, X2)dx = 0 , 

B 2 


if L + P < M. We obtain 




□ 


Proof (of Theorem 3.3): We use the same notation as in the proof of Theorem 3.2 for the functions Ui = 
(ail) ^ W-i = u and the discontinuity curves 7 t, t = — 1,1 and the outer normal n''(t) of 7 ^ at 

t for i = 1, —L We decompose our index set A. The sets A° and A^, Af are defined exactly as in the proof 
of Theorem 3.2 We do, however, decompose A^ differently, into A^’“, A^’^ and A^’'^, to be defined below. 

We denote by A^’“ the indices in A^ := A \ (A° U A^), such that for (j, 2^/^s , m, t) G A there exists some 
t G supp '0j\22/2s,m,(,fl7t such that we have G (tanpo —3-2“^/^,tanpo+3-2“^/^) and n‘'{t) = ±{cos po,sm po) 
for some po G (—7r/2, 71 / 2 ). 

Additionally, there are non-aligned shearlet elements, where the angle between the shearlet and the 
discontinuity curve is less than 7 r /4 and that touch the discontinuity curve, we will call these indices A^’^. 
More precisely, these are indices (j, 2'1/^s, m, t) G A that are not in A° U A^ U A^’“ and for which there 
exists some t G supp ipj^ 2 ^/ 2 ^ m,i h. such that we have that (i[- 7 r/ 2 ,- 7 r/ 2 ](arctan(s'’), po) < ’’"/d and n‘'(t) = 
±(cospo,sinpo) for some po G (— 7 r/ 2 , 7 r/ 2 ). 

Lastly, there are indices, such that the angle between the curve and the shearlet is larger than 7 r/ 4 , i.e. 
(j, 2 '^/^s,to, t) G a \ (A° U A^ U A^’“ U A^’^) such that d[_,n._.n.](arctan(s''), po) > '^/4 we will call these indices 

For j G Z we denote by Aj, A^’“, A^’**, A^’^ the indices of the respective index set, with scale equal to j. 
We have that #{(j, A:,m, t) G Aj : (u, ■!/’pfc,m,t) 0} ^ 2^-^. Furthermore, 


|A'’“|<2J, \Af\<2\ |A^'’"|<2iL 

Again we consider three different parts: 

Part 1: A°: 

Since |A°| < 00 we certainly have 

n>N 

where c(it)* denotes a non-increasing rearrangement of (|(u,' 0 j,fe,m,t)P)(j,fe,m,t)GAo. 







Part 2: supp ’tpj,k,m,L n Ut={-i 4 } = ^ ■ 

Using we obtain that 

I (^7 I ^2 I ('a^, |- 

Since G C°° and 9 has M > I + 9/2 vanishing moments in xi-direction and supp 9j^k,m,L is of length 2~^ 
in the direction indicated by b we can estimate 

< 2 - 2 ('+l)U 

Let e > 0 and p = {I + 9/4 — e)~^. From the cardinality of Aj we have that 

o, (‘+ 1 ) 

Ai jGN 


Using Stechkin’s lemma Q, we obtain that X]n>Arc('u)* < N 2(i+7/4 e)^ where c{u)* denotes a non¬ 
increasing rearrangement of (|(M,'!/'j,fe.m.JP)(i,fc,m,OGAi- 
Part 3: supp n Ut={-i 4 } ^ 0 : 

As already established in the proof of Theorem 3.2 we have for (j, k, m, t) € A^’“ 
and for {j, k, m, l) G A^’^ we have 

We continue with the coefficients of First of all we invoke (§ and obtain 

I (^j I ^2 ^ 9jk^m,i,)\- 


Next we want to estimate |(ut, dyfe,m,t)l using Lemma 3.4 Since G £°°{v) we have that 

Uu= gi+ Xb92, with 31,52 G 


for all P G N and in particular for P = \l + . In addition we have that dD G C° °. Furthermore 9 G 

has M >2L vanishing moments and L > \l + ^~\. Hence we can apply Lemma 


3.4 


to obtain that 


(yj^k,771,1/ 




Consequently 


< 2 -'^ 2 -('+T)i = 2 i('+t). 


Invoking the sizes of A^’“, A^’^, A^’'^ we can compute, that 


I {u, V’i.fc.m, J r < XI I + X I (“> T + X I T 

A2 A2.“ A^.t AS,': 

< X] 2'^2“^^2“'"^^-^ -h X] 25 2“^*--I- X] 2^2“5P(^+lb < oo. 

jeN jgn i'gn 


(9) 


Another application of Stechkin’s lemma yields that X]n>Ar *^W)n < Xf 2(i+7/4 e)^ where c(u)* denotes a 
non-increasing rearrangement of {\{u,ipj^k,m,/)\‘^){j,k,m,L.)GA^■ Combining the estimates from Part 1, 2, and 3 
yields with the dual frame characterization 0 that 

Ik- X k>V’A)V’AllL(RS) ^ X + 

where c(4t)* denotes a non-increasing rearrangement of (|(m, V'a)P)agA- □ 
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